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Abstract 

Icing scaling tests were performed in the NASA Glenn Icing Research Tunnel to evaluate a new 
scaling method (Ref. 1), developed and proposed by Feo for glaze icing, in which the scale liquid water 
content and velocity were found by matching reference and scale values of the nondimensional water-film 
thickness expression and the film Weber number. For comparison purpose, tests were also conducted 
using the constant We L method for velocity scaling. The reference tests used a full-span, fiberglass, 
91.4-cm-chord NACA 0012 model with velocities of76and 1 00 knot and MVD sizes of 150 and 195 pm. 
Scale-to-reference model size ratio was 1 :2.6. All tests were made at 0° AOA. Results will be presented 
for stagnation point freezing fractions of 0.3 and 0.5. 

Introduction 

The ability to conduct icing tests with subscale aircraft components in an icing wind tunnel depends 
critically on our understanding of physical mechanisms that affect the ice accretion process and our 
ability to measure and control test conditions in an icing tunnel. The need to develop appropriate scaling 
methods has been long recognized, and research efforts have begun as early as in the 1950s and continue 
to the present. Recently, Anderson (Ref. 2) and Anderson and Tsao (Ref. 3) have completed two NASA 
reports that provide a detailed technical review of recommended scaling methods for ice accretion on 
unprotected, unswept aerodynamic surfaces in Appendix C and super-cooled large droplets (i.e., SFD) 
conditions, respectively. 

It was concluded from those two references that acceptable scaling results could be achieved by 
matching the p 0 > A, n 0 , and We/ . With scale model size selected, by matching scale and reference values 
of We/, the scale velocity can be determined. By matching (3 0 the scale MVD can be found. Reference 2 
also showed that the effects of temperature and L WC are not independent, but interact through the 
freezing fraction. Therefore, with scale LWC chosen, by matching n 0 the scale temperature can be 
calculated. Finally, by matching A c the scale time can be established. For the scale test, then, only 
temperature, velocity, MVD and time have to be calculated from the known (reference) values of the 
similarity parameters. 

The results presented here are part of an effort to develop scaling methods for super-cooled large 
droplets in glaze icing conditions. During glaze-ice accretion water does not freeze immediately in the 
zone of impingement. The characteristics of the resulting layer of liquid surface water are believed to 
influence the shape of the accreted ice to some extent and are the subject of a number of ongoing studies. 
Bilanin (Ref. 4) and Bilanin and Anderson (Ref. 5) advocated adding another similarity parameter, the 
Weber number, to scaling requirements to address surface-water effects. Two forms of the Weber number 
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that might be applied include that based on the droplet MVD and water properties (We) and one based on 
model size and air properties (We c ). A third Weber number was suggested by Kind (Ref. 6); this one was 
based on the thickness of the water film at the leading edge and air properties (We,). 

Anderson and Ruff (Ref 7) and Anderson (Ref 8) reported the best scaling results for glaze ice were 
achieved when a compromise between constant We and constant Re was used with the modified Ruff 
method. A possible explanation for this result comes from studies of droplet-surface impact (Ref. 9) for 
no icing applications. These studies correlated splashing effects with a K factor that depended only on We 
and Re s . Other surface-tension effects such as water-film breakup and rivulet formation are also We and 
Re dependent. 

Finally, in a recent study at above-freezing temperatures, Feo and Urdiales (Ref. 10) measured the 
water-film thickness, h, for heavy-rain conditions. The nondimensional film thickness, h/d, was found to 
correlate with We and Re 5 . The authors suggested that h/d might itself serve as an appropriate similarity 
parameter for scaling studies. Although spray characteristics for that study were very different from 
Appendix-C conditions, the heavy-rain h/d was evaluated as a similarity parameter in studies comparing 
how methods to select scale velocity affected ice shape (Refs. 7 and 8). However, matching the scale and 
reference h/d of Feo and Urdiales gave scale ice shapes that were generally a poorer match of the 
reference than did other methods of finding scale velocity. Later, Feo (Ref. 11) measured the thickness of 
the water film for MVD and LWC conditions near those in an icing cloud and correlated the results with 
LWC and Re. In the study he suggested another Weber number based on the water film thickness and 
water properties (We/,) as a similarity parameter for scaling studies. Evaluations of the constant-water- 
film-thickness and constant-Weber-number methods in scaling tests using this new h/d correlation have 
been made (Ref. 12) and the results compared reasonably well with those using the Weber number We. 

In another more recent study (Ref. 1), similarity parameters involved in surface water film dynamics 
were proposed. The film Weber Number (We/) that was defined before (Ref. 13) was introduced with the 
nondimensional film thickness, being both required to be matched at scale and reference conditions. The 
parametric ( h/d) expressions were obtained from two sets of experiments that were conducted in one of 
INTA’S Low Speed Tunnels; one being for Appendix C (Refs. 1 1 and 14) and the other for SLD droplets 
(Ref. 15). In addition, Feo has also suggested a slightly different approach from the modified Ruff method 
as recommended in References 2 and 3 to determine appropriate scale test conditions. Three of the five 
similarity expressions used were from well-established methods (droplets’ trajectories, freezing fraction 
and ice accumulation) and the other two requirements were obtained by using We/= const, and (h/d) = 
const., the last one being taken from experiments. Depending on the expression of (h/d) used; we have 
results either for Appendix C or SLD droplets. However there was no exact icing test results available at 
the time to properly evaluate the proposed scaling method and the need of more experimental work 
related to this specific method was recommended in Reference 1 . 

This evaluation was a goal of the present study, along with the evaluation of the constant Wei method 
in glaze icing conditions with SLD reference droplet sizes in the range of 1 00 to 200 pm. 

Nomenclature 

A c Accumulation parameter, dimensionless 
b Relative heat factor, dimensionless 

c Airfoil chord, cm 

d Cylinder radius or twice the leading-edge radius of airfoil, cm 

h Water film thickness, cm 

K 0 Modified inertia parameter, dimensionless 

L WC Cloud liquid-water content, g/m 3 

MVD Water droplet median volume diameter, pm 

n Local freezing fraction, dimensionless 

77o Stagnation point freezing fraction, dimensionless 

Re Reynolds number of model, dimensionless 
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Reg Reynolds number of water drop, dimensionless 
t st Static temperature, °C 

V Air velocity, kt 

We Weber number based on droplet size and water properties, dimensionless 
We c Weber number based on model size and air properties, dimensionless 

We/, Weber number based on water- film thickness, air velocity and water properties, dimensionless 

We/ Weber number based on water- film thickness, water properties and film velocity, dimensionless 

Wei Weber number based on model size and water properties, dimensionless 

We, Weber number based on water-film thickness and air velocity and properties, dimensionless 

Weg Weber number based on droplet size and air properties, dimensionless 

[to Collection efficiency at stagnation line, dimensionless 

<|) Droplet energy transfer parameter, °C 

x Accretion time, min 

Subscripts 

R reference 
S scale 
/ water film 
st static 
tot total 


Similarity Parameters 

The similarity parameters used in this study were based on the work originally done by Ruff 
(Ref. 16). The current scaling method (i.e., the modified Ruff method) required matching scale and 
reference values of the key similarity parameters, Po, A c , no, and Wei with scale LWC chosen at user’s 
convenience. In Feo’s approach, matching reference and scale film Weber number Wc, sct the scale 
velocity, and matching nondimensional film thickness parameter h/d set the scale LWC. The equations for 
the similarity parameters will be presented here without much discussion. Therefore, readers who are 
interested in the physical descriptions and detailed derivations of these parameters are referred to 
References 1 to 3 and the references given therein. 

To maintain the droplet trajectory similitude, Langmuir and Blodgett (Ref. 17) introduced the 
modified inertia parameter, K 0 , defined as 



L Stokes 


v 8 / 


(1) 


In Equation (1), AT is the inertia parameter, 


P W MVD 2 V 
18 d p 


( 2 ) 


Departing slightly from Langmuir and Blodgett in this study, d is twice the leading-edge radius of 
curvature for airfoils. For the NACA 0012 airfoil model, a leading-edge radius of 0.0158c was used (see 
Abbott and von Doenhoff (Ref. 18)), where c is the airfoil chord. L/L St0 kes is the droplet range parameter, 
defined as the ratio of actual droplet range to that if Stokes drag law for solid-spheres applied. It is a 
function only of the droplet Reynolds number, Re 6 . 

VMVDp 
Res = 
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This study used a curve fit to Langmuir and Blodgett’s tabulation of the range parameter as given in 
the following expression: 


1 1 

^Stokes (0.8388 + 0.001483 Re § +0.18477^") 


( 4 ) 


Of more practical interest than K 0 is the collection efficiency at the stagnation point, p 0 , which was 
shown by Langmuir and Blodgett to be a function only of K 0 , 


Po = - 


1.401 K 0 - 


.84 


1+1.401 


.84 


(5) 


Thus the droplet trajectory similarity is satisfied if K 0 , s = K 0 , R (so does Po,s = Po,«), and the scale drop 
size, i.e., scale MVD, is determined. To ensure water-catch similarity, the accumulation parameter is 
introduced: 


LWCVx 

dpi 


(6) 


If all the water impinging on the leading edge freezes at that location and the leading-edge collection 
efficiency is 100 percent, A c directly becomes a measure of the normalized thickness of ice that will 
accrete. The scale accretion time can be found from A c , s = A C , R . However if it is not possible to find scaled 
conditions that permit a match of K 0 (and therefore p 0 ) and A c separately, it is recommended that the 
product of P oA c be matched provided the two collection efficiencies, p 0 ,.s and p 0 , fo be within 10 percent, 
see Reference 2 for details. 

The rate at which the water freezes on a surface depends on the magnitude of local heat transfer 
imbalance. For glaze ice, it is known that the fraction of water that freezes is less than unity, and the 
motion of unfrozen surface water can have an effect on the resulting ice shape. Therefore, it is important 
to maintain surface energy and surface-water dynamics similarities for glaze ice accretions. The freezing 
fraction is formally defined as the ratio of the amount of water that freezes at a given surface location to 
the total amount of water that impinges at that location. From Messinger’s (Ref. 19) steady-state surface 
energy balance analysis, the stagnation point freezing fraction can be written as 


«o = 


Cp,WS 




O' 

by 


(7) 


The key terms in this formulation include (j) and 9, which have dimensions of temperature and relate 
to the water drop energy transfer and air energy transfer, and b, the relative heat factor, which was first 
introduced by Tribus, et al. (Ref. 20) 


9 t f t st 


V 2 


( 8 ) 
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LWCV?, 0 c P ' WS 

K 


( 10 ) 


Equation (9) from Ruff (Ref. 16) has included compressibility effects. Various incompressible forms 
of 0 have also been used by Charpin and Fasso (Ref. 21) and others; however, the differences are not 
significant mainly due to the fact that, for most icing conditions, the Mach number is relatively low. Since 
only limited combinations of MVDs and LWCs in the SLD regime are currently available in the IRT, the 
stagnation point freezing fraction is matched to find the scale static temperature instead of the scale L WC 
value in this study. 

Because the original Ruff method does not restrict the value of scale velocity, an additional similarity 
parameter can be used to determine Vs. In 2003 Anderson and Tsao (Ref. 22) had provided experimental 
evidence from past studies to show that a similarity parameter dependent on the ratio V x c y /a : must be 
included in scaling methodology to account for surface-water dynamics effect in glaze ice accretions, 
although the powers x, y, and z are not yet determined. The length may not be chord itself but rather some 
physical characteristic L related to chord; for example, the water-film thickness. Likewise, the velocity 
could also be of the water-film which is related to V. Thus a Weber number based on L and V 


We L = 


v 2 l Pw 


( 11 ) 


has been suggested as a potential additional similarity parameter to supplement Ruffs basic scaling 
method. Studies by Bartlett (Refs. 23 and 24) and Oleskiw, et al., (Ref. 25) found no measurable effect of 
pressure on ice shape. These observations suggest that water density is a better choice than air density for 
Equation (11). In this study the Wei is based on the twice the nose radius of the airfoil: 


We z = 


V 2 d p w 


( 12 ) 


with the understanding that L cc d. The scale velocity found from matching Wei, 5 = Wc/,/ ( is 


v s =v R 


r d S V/2 


\ d s ) 


(13) 


Another approach presented by Feo (Ref. 13) involved matching simultaneously the water-film 
thickness and a Weber number that used the water-film velocity and water-film thickness: 


We f = 


v f h fPy 


( 14 ) 


By equating shear stresses at the air-water interface, the water-film velocity can be written as 
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V f = Re 


1/2 




V d J 


V 


(15) 


When Equation (15) is substituted into Equation (14) and letting scale and reference values of (h f /d) 
be equated while matching Wc/, iS - = Wc /5ff , the expression becomes 


(ReWeJ 5 =(ReWe^ 
From Equation (16) the scale velocity for this Feo method is 


(16) 


v s = v R 


f d ^ 2/3 

a R 

\ d S ) 


(17) 


The Reynolds number Re shown above is also based on the twice the nose radius of the airfoil: 


Re = 


Vd p 


T 


(18) 


For large SLD drop sizes where droplet impact could strongly interact with the surface air-driven film 
flow, the best correlation of the experimental film thickness data (Ref. 1 5) was found to be 


hid <x 


f \l/2 

' LWC' 


V Fw J 


Re- V4 We c - 1/2 (We 5 1/2 Re 5 - 1 ) 


(19) 


By matching scale and reference values of nondimensional film thickness ( h/d '), a seal eLWC was 
determined 


LWC s =LWC r 


Ug_ 

\ d s j 


( 20 ) 


Test Description 

Facility, Model, and Procedures 

The scaling tests were performed in the NASA Glenn Icing Research Tunnel (1RT). The 1RT is a 
closed-loop, refrigerated, sea level tunnel with a 1.8 by 2.7 m rectangular test section. The icing cloud is 
generated by operating 10 spray bars, a configuration in use since 1998. 

The 1RT cloud calibrations for both Appendix-C and SLD conditions used for these tests were 
performed in the summer of 2008. The LWC measurements were made using icing blade method as 
reported previously (Ref. 26). Flowever, during the 2008 calibration there was not enough time to perform 
MVD measurements. The MVDs reported in this paper are based on an analysis of the MVD calibration 
data completed in February 2006. In addition, because only a few specific MVD-LWC combinations at 
speeds of 100, 150, 200, and 250 kt (i.e., 51, 77, 103, and 128 m/s) have been calibrated to date in the 
SLD regime, additional LWC measurements for SLD conditions were made in the 1RT on September 22, 
2008, with an icing blade. Therefore SLD tests are constrained to these particular conditions. 
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The models used were NACA 0012 airfoil sections with chords of 91.4 and 35.6 cm. The 91.4-cm- 
chord airfoil is pictured in Figure 1(a). It was a full-span, fiberglass model at 0° angle of attack and served 
as the reference model. The 35.6-cm-chord scale model was of 61 -cm span and made of aluminum. It was 
mounted vertically between splitter plates at the center of the 1RT test section as shown in Figure 1(b). 
Horizontal lines at the leading edge were drawn at the tunnel vertical center (model mid span) and 
±2.5 cm from the center to locate ice-tracing templates. Vertical lines were also placed at increments of 
2.5 cm (labeled in inches on the model), measured along the surface from the stagnation line. These 
marks helped to identify sites on the model for close-up photographs of feather structure details. Because 
of the quick-start capability of the 1RT spray system, no shielding of the models was required during the 
initiation of the spray. 

In preparing for a test, the temperature and airspeed in the test section and the air and water pressures 
on the spray manifolds were set. When these conditions had stabilized, the spray nozzle valves were 
opened to initiate the spray. The spray was timed for the required duration, and then turned off. The fan 
was brought to a full stop and the researchers entered the test section to document the ice shape with hand 
tracings. Close-up photographs were also taken with a hand-held digital camera. 

To record the ice shapes, a thin slice was first melted through the ice normal to the model surface. A 
cardboard template was then placed into this slit and an outline of the ice shape traced by pencil, giving a 
two-dimensional cross section of the ice. Tracings were taken at the vertical center of the tunnel (91 cm 
from the floor) and at 2.5 cm above the center. The ice shapes so recorded were digitized using an 
automated line-following feature in the image-analysis software, SigmaScan Pro (Ref. 27). The results 
presented in this study are from 1RT test entry in September 2008. Since the shape differences between 
the two tracing locations were never significant, only centerline shapes will be reported here. 



(a) 91.4-cm-Chord NACA 0012 Model (b) 35.6-cm-Chord NACA 0012 Model 


Figure 1. — Models installed in IRT test section. 
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Uncertainty Analysis 

Estimates of the uncertainty in the reported average conditions were made by considering inherent 
errors of instruments, temporal fluctuation and spatial variation of the instrument readings in the test 
section, and uncertainty in tunnel calibration of MVD and L WC. Recorded air temperature was believed to 
be accurate to ±0.2 °C, although variations during the period of an icing spray increases the uncertainty 
for reported average temperatures to about ±0.5 °C. The uncertainty in air velocity was estimated to be 
±2 kt. For Appendix-C conditions the net uncertainty in MVD was estimated at ±12 percent. For SLD 
conditions it may have been as much as ±20 percent. These uncertainties are not referenced to an absolute 
value of MVD, which is u nk nown. Repeatability and scatter in the L WC calibration data suggests the 
uncertainty is about ±12 percent for both Appendix-C and SLD conditions. 

The test-parameter uncertainties were used to estimate the following uncertainties in the similarity 
parameters for the SLD tests the uncertainties were: 2 percent in p 0 , 12 percent in A c , 10 percent in « 0 , 

2 percent in Re, 3 percent in Wei, and 5 percent in We/. 

Results 

In planning the test matrix, reference test conditions were chosen so that both reference and scale 
conditions would fall within the 1RT operating envelope. Due to very limited 1RT test time available for 
the scaling evaluation this study is only focused on the SLD droplets. Also, since only limited 
combinations of MFDs and LWC s in the SLD regime are currently available in the 1RT, it was decided 
early on for time-saving purpose that the scale LWC values obtained from Equation (20), as recommended 
by Feo to ensure constant h/d, would also be used for the evaluation of the constant Wei method. 

Figures 2 to 9 present the nondimensional reference and scale ice shapes recorded using two methods 
for finding scale velocity: the constant Wei and constant We/ for stagnation point freezing fraction of 0.3 
and 0.5 at reference velocities of 76 and 100 kt. For each figure, reference ice shape was shown shaded, 
while a solid line indicated the scale shape. The table below each figure gave the test conditions and 
similarity parameters for each pair of reference and scale tests. The conditions given were the average 
conditions recorded over the duration of each test, which can sometimes differ slightly from the planned 
set points. The parameters in the tables were calculated from these average conditions. 

Stagnation Point Freezing Fraction of 0.3 

Figures 2 and 3 show reference and scale ice shape comparisons for the constant We L and constant We/ 
methods, respectively. The same reference ice shape was used for both comparisons. The reference model size, 
velocity and MVD were 91.4 cm, 76 kt, and 150 pm. The scale model size was 35.6 cm. The coordinates of the 
ice shapes were all normalized by the corresponding model chord. 

In Figure 2 the scale and reference values of Po, »o, and Wei matched within 3 percent and A, was just 
within 9 percent. This was because the reference LWC value was reduced to 87 percent of the planned value 
(i.e., 1.68 instead of 1.93 g/m 3 ) and the scale LWC was reduced by 5 percent from its planned value (i.e., 0.95 
instead of 1.00 g/m 3 ) based on the calibration performed in the 1RT on September 22, 2008. The scale ice 
shape however closely simulated the reference main ice shape and feather region, even including smaller 
feathers further aft on the surface. 

In Figure 3 the scale and reference values of Po matched well but A c and Po A c were just within 1 6 percent, n 0 
was within 10 percent, We/agreed only within about 37 percent and h/d was within 15 percent. Such large 
disagreement of We/and h/d occurred mainly because, in addition to the decrease in the reference LWC value 
mentioned previously, now the scale L WC value was further increased to 1 17 percent of the planned value (i.e., 
1.17 instead of 1 .00 g/m 3 ) based on the calibration data. The Olsen method was used in this case to adjust the 
spray time and air total temperature for scale test due to the LWC value change from the calibration such that « 0 
and Po A c would stay close to the planned value. As a result, the size and shape of the reference ice were simulated 
fairly well by the scale test in the leading-edge region, though missing the hom angle in the upper feather region. 
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Figure 2. — Scaling from 91.4 to 35.6-cm-chord with We/, matched. NACA 0012 airfoils; n 0 , 0.3; V R , 76 kt. 
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Figure 3. — Scaling from 91.4 to 35.6-cm-chord with We f and ( h/d ) matched. NACA 0012 airfoils; n 0 , 0.3; V R , 76 kt. 
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In addition, Figures 4 and 5 gave SLD scaling results for a higher reference velocity of 100 kt with 
the 91.4-cm-chord 195-pm-MFD reference conditions scaled to 35.6 cm. In both cases, the scale and 
reference values of Po, Po A c and n 0 matched within 3 percent. The reference and scale Wei in Figure 4 
matched well and the We/ and h/d in Figure 5 agreed within about 3 percent. The 35.6-cm-chord scale 
tests produced main ice shapes in fairly good agreement with the reference though the constant We L 
method did seem to simulate the leading-edge region better. Even the sizes of large feathers adjacent to 
the main shape and the smaller feathers further aft were simulated reasonably well (except for a few 
larger feathers further aft on the lower surface in both figures). 



c = 91.4 cm, 09-26-08 Run 1 
c = 35.6 cm, 09-30-08 Run 2 


Date/Run 

C, 

hob 

V, 

MVD , 

LWC , 

T, 

Po, 

A c 

Pdc 

”o 

We L , 

We, 

h/d. 


cm 

°c 

kt 

jam 

g/m 3 

min 

percent 




10 6 

1 0 -16 

io- 9 

09-26-08/01 

91.4 

-7.8 

100 

198 

1.08 

14.0 

96.3 

1.75 

1.69 

0.31 

1.17 

0.74 

0.82 

09-30-08/02 

35.6 

-2.9 

159 

87 

0.57 

6.62 

96.2 

1.80 

1.73 

0.31 

1.16 

0.86 

1.04 


Figure 4. — Scaling from 91.4 to 35.6-cm-chord with Wei matched. NACA 0012 airfoils; no, 0.3; Vr, 100 kt. 



c = 91.4 cm, 09-26-08 Run 1 
c = 35.6 cm, 09-30-08 Run 1 


Date/Run 

C, 

hob 

V . , 

MVD, 

LWC, 

T, 

Po, 

A c 

PckT 

”0 

Wei, 

We, 

h/d. 


cm 

°c 

kt 

Jim 

g/m 3 

min 

percent 




10 6 

10” 6 

10“ 9 

09-26-08/01 

91.4 

-7.8 

100 

198 

1.08 

14.0 

96.3 

1.75 

1.69 

0.31 

1.17 

0.74 

0.82 

09-30-08/01 

35.6 

-2.7 

185 

86 

0.54 

5.65 

96.4 

1.69 

1.63 

0.32 

1.56 

0.72 

0.84 


Figure 5. — Scaling from 91.4 to 35.6-cm-chord with We^ and (h/d) matched. NACA 0012 airfoils; no, 0.3; Vr, 100 kt. 
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Stagnation Point Freezing Fraction of 0.5 

Figures 6 and 7 each compare reference and scale ice shapes with the 91.4-cm reference conditions 
scaled to 35.6 cm. The same reference ice shape was used for both comparisons. The reference velocity 
and MVD were 76 kt and 150 pm. In Figure 6 the scale velocity was calculated from Equation (13), 
whereas that for Figure 7 was obtained from Equation (17). 



x/c 


c = 91.4 cm, 09-26-08 Run 4 
c = 35.6 cm, 09-30-08 Run 8 


Date/Run 

C, 

cm 

hot > 

°c 

V, 

kt 

MVD , 
(dm 

LWC, 

g/m 3 

T, 

min 

Po, 

percent 

A c 

P (Ac 

»o 

W ei , 

10 6 

We/, 
10 36 

h/d, 

io - 9 

09-26-08/04 

91.4 

-18.1 

76 

148 

1.78 

11.2 

94.1 

1.76 

1.66 

0.52 

0.68 

3.03 

1.68 

09-30-08/08 

35.6 

-9.0 

121 

68 

0.94 

5.26 

94.2 

1.79 

1.69 

0.52 

0.67 

3.29 

2.09 


Figure 6. — Scaling from 91.4 to 35.6-cm-chord with Wei. matched. NACA 0012 airfoils; n 0 , 0.5; V R , 76 kt. 



x/c 


c = 91.4 cm, 09-26-08 Run 4 
c = 35.6 cm, 09-30-08 Run 9 


Date/Run 

C, 

hot •> 

V, 

MVD , 

LWC, x, 

Po, 

A c 

P oA c 

»o 

Wei, 

We,, 

h/d. 


cm 

°c 

kt 

(im 

g/m 3 min 

percent 




10 6 

10“ 36 

IO " 9 

09-26-08/04 

91.4 

-18.1 

76 

148 

1.78 

11.2 

94.1 

1.76 

1.66 

0.52 

0.68 

3.03 

1.68 

09-30-08/09 

35.6 

-10.4 

141 

66 

1.18 

4.00 

94.3 

1.99 

1.88 

0.48 

0.91 

4.37 

1.97 


Figure 7. — Scaling from 91.4 to 35.6-cm-chord with Wef and ( h/d ) matched. NACA 0012 airfoils; no, 0.5; Vr, 76 kt. 
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In Figure 6 the scale and reference values of p 0 , Po 4„ n 0 and We L all matched within 2 percent, and 
the size and shape of the reference ice were well simulated by the scale test in the leading-edge region. 
However, the scale ice feathers showed relatively larger formations than recorded in the reference test. 
Further study to isolate these large-feather conditions is needed to improve scaling performance. In 
Figure 7 the scale and reference values of p 0 matched well but A, and P oA c were just within 12 percent, n 0 
was within 8 percent, We/ agreed only within about 31 percent and h/d was within 15 percent. The large 
differences in matching We/ and h/d occurred because the scale LWC value was increased to 1 18 percent 
of the planned value (i.e., 1.18 instead of 1.00 g/m 3 ) while the reference LWC value was reduced by 8 
percent from its planned value (i.e., 1.78 instead of 1.94 g/m 3 ) based on the calibration data. The Olsen 
method was again applied to adjust the spray time and air temperature for scale test to keep n 0 and P oA c 
close to the planned value. The shape of the reference ice was simulated fairly well by the scale test in the 
leading-edge region, but the size of the scale ice including the horns was bigger than the reference. In 
addition, many large ice feathers were observed for scale test in Figure 7. 

For a higher reference velocity of 1 00 kt, similar ice shape comparisons were made with the 9 1 ,4-cm 
reference conditions scaled to 35.6 cm, and the results were shown in Figures 8 and 9. The reference 
MVD was 195 pm and the scale MVD was about 87 pm. In both cases, the scale and reference values of 
Po, Po A c , and n 0 matched within 3 percent. Reference and scale We/ in Figure 8 matched within about 
2 percent, and the We/ and h/d for Figure 9 agreed within 3 percent. In Figure 8 the scale ice shape was 
able to simulate the reference main ice shape and feather region well. In Figure 9, however, the scale test 
produced a smaller horn angle than the reference, although it had the correct horn size and reproduced 
well the feather size and density. 

Given the repeatability variations encountered in the IRT for SLD conditions, it can be concluded that 
for SLD droplets in conditions where the droplet impact could interact strongly with the surface shear- 
driven film flow the constant (We L and h/d) method seems to be a better choice of similarity parameters 
for determining scale velocity and liquid water content than the constant (We/ and h/d) method at 
stagnation point freezing fractions of 0.3 and 0.5. Additional tests at these freezing fractions are needed to 
better assess the constant (We^ and h/d) method. The preliminary assessment is, however, that any of the 
two methods can produce acceptable scaling simulation for « 0 = 0.3 and 0.5. If scale simulations are 
needed at such low freezing fractions, it is recommended that, due to the variability of shapes, tests be 
repeated with more than one test entry. 



x/c 


c = 91.4 cm, 09-26-08 Run 3 
c = 35.6 cm, 09-30-08 Run 5 


Date/Run 

C, 

cm 

hot > 

°c 

V, 

kt 

MVD , 
jam 

LWC, 

g/m 3 

X, 

min 

Po, 

percent 

A c 

P (Ac 

no 

W ei , 

10 6 

We r , 
10 36 

h/d, 

10~ 9 

09-26-08/03 

91.4 

-14.1 

100 

198 

1.07 

14.0 

96.3 

1.75 

1.69 

0.52 

1.18 

0.74 

0.81 

09-30-08/05 

35.6 

-5.9 

159 

87 

0.57 

6.62 

96.2 

1.80 

1.73 

0.51 

1.16 

0.87 

1.03 


Figure 8. — Scaling from 91.4 to 35.6-cm-chord with WeL matched. NACA 0012 airfoils; no, 0.5; Vr, 100 kt. 
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x/c 


c = 91.4 cm, 09-26-08 Run 3 
c = 35.6 cm, 09-30-08 Run 6 


Date/Run 

C, 

t toti 

V, 

MVD , 

LWC, 

T, 

Po, 

A c 

Pc Ac 

«o 

Wei, 

We,, 

h/d. 


cm 

°c 

kt 

jam 

g/m 3 

min 

percent 




10 6 

10‘‘ 6 

10' 9 

09-26-08/03 

91.4 

-14.1 

100 

198 

1.07 

14.0 

96.3 

1.75 

1.69 

0.52 

1.18 

0.74 

0.81 

09-30-08/06 

35.6 

-6.0 

186 

86 

0.54 

5.65 

96.4 

1.70 

1.64 

0.53 

1.58 

0.72 

0.83 


Figure 9. — Scaling from 91.4 to 35.6-cm-chord with Wet and (h/d) matched. NACA 0012 airfoils; no, 0.5; Vr, 100 kt. 


Conclusion 

Icing scaling tests were performed in the NASA Glenn Icing Research Tunnel to evaluate a scaling 
method, developed and proposed by Feo for glaze icing conditions, in which the scale liquid water 
content and velocity were found by matching reference and scale values of the INTA nondimensional 
water- film thickness expression h/d and the film Weber number We/. For comparison, tests were also 
performed using the constant Wei method. The reference tests used a lull-span, fiberglass, 91.4-cm-chord 
NACA 0012 model with velocities of 76 and 100 knot and MVD sizes of 150 and 195 pm. Scale-to- 
reference model size ratio was 1:2.6. All tests were made at 0° AOA. Results were presented for 
stagnation point freezing fractions of 0.3 and 0.5. 

For freezing fraction of 0.5, all methods resulted in scale accretions closely matching the size of the 
reference, but somewhat better agreement with the reference hom angles was typically achieved with the 
constant (Wei and h/d) method than with the constant (We ; and h/d) method. At a freezing fraction of 0.3, 
the two methods were judged equal in providing simulations of the reference shape within the ice-shape 
repeatability. At low freezing fractions, poorer repeatability has been observed than for higher freezing 
fractions. Conclusions of this study regarding the constant-Wei and constant-We/ methods of finding 
scale velocity agreed with those of Reference 12 and, in which any scaling method considered for glaze 
ice accretion has to be consistent with the presumption that both Re and We are important parameters in 
describing those water-film phenomena that affect ice accretion. 
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